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Abstract 

This study investigates Parker instability in an interstellar medium (ISM) near the Galactic plane using three- 
dimensional magneto-hydrodynamic simulations. Parker instability arises from the presence of a magnetic field in 
a plasma, wherein the magnetic buoyant pressure expels the gas and cause the gas to move along the field lines. The 
process is thought to induce the formation of giant molecular clouds in the Galaxy. In this study, the eff'ects of cosmic- 
ray (CR) diff'usion are examined. The ISM at equilibrium is assumed to comprise a plasma fluid and a CR fluid at 
various temperatures, with a uniform magnetic field passing through it in the azimuthal direction of the Galactic disk. 
After a small perturbation, the unstable gas aggregates at the footpoint of the magnetic fields and forms dense blobs. 
The growth rate of the instability increases with the strength of the CR diff'usion. The formation of dense clouds 
is enhanced by the eff'ect of cosmic rays (CRs), whereas the shape of the clouds depends sensitively on the initial 
conditions of perturbation. 
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1. Introduction 

The gaseous disk of our Galaxy has been suggested to 
be unstable in response to perturbations in the magnetic 
field line that lies parallel to the Galactic plane. The dy- 
namic process of the subsequent mixing of gas is called 
Parker instability in the astronomical literature (Parker 
1966) [1]. The morphology of this instability can ex- 
hibit two modes - undular mode and interchange mode 
(Kim et al. 1998) [2]. The undular mode propagates 
along the magnetic field in the x direction and is capa- 
ble of exciting the gas in the vertical direction. An ini- 
tially small perturbation in the undular mode can distort 
the uniform magnetic field and cause convection at full 
scales (Lee 2009) [3]. However, the interchange mode 
travels in the y direction, and may or may not aff'ect the 
horizontal field lines. Hence, when the magnetic field 
has an oscillating component in the Galactic disk plane, 
convective instability may occur. 

Once the gas rises because of magnetic buoyancy, the 
plasma can only fall along the magnetic field lines if 
the magnetic field is frozen into the plasma. Thus, a 
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magnetic arch that protrudes from the bottom horizon- 
tal magnetic field is associated with dense footpoints. 
The key observation that supports Parker instability is 
that of giant molecular loops close to the Galactic cen- 
ter (Fukui et al. 2006) [4]. The gas aggregated at the 
footpoint is expected to collapse eventually, becoming 
a giant molecular cloud and a region that hosts stellar 
formations. The undular mode is also thought to power 
the magnetic dynamo, which sustain the magnetic field 
in the Galaxy (Hanasz et al. 2004) [5]. 

In this work, CR diffusion is adopted to elucidate 
Parker instability. Cosmic rays are highly energetic par- 
ticles, 90% of which are protons. The speed of indi- 
vidual CRs particles is close to that of light, but the 
bulk motion of CR is diff'usive and the bulk speed of 
CRs is of the order of the Alfven speed. The CR parti- 
cles propagate along the magnetic field lines as long as 
their gyroradius is significantly smaller than the charac- 
teristic spatial scales of the magnetic field. Cosmic rays 
may critically aff'ect the dynamics of the ISM, since the 
energy density of CRs is of the same order as that of 
the magnetic field and turbulent gas motions. The im- 
portance of the eff'ects of CRs has been acknowledged, 
and various works on CRs and Parker instability have 
been published. Hanasz et al (2004) identified the bene- 
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fit of considering the CRs eff'ect in the Galactic dynamo 
problem. Otmianowska-Mazur et al. (2009) [6] con- 
structed a numerical model for the CR driven dynamo; 
they showed that excitement of fast growth in the mag- 
netic flux and energy can account for the X- shaped mag- 
netic structures of the galactic halos in several galaxies 
(Beck 2009) [7], which cannot be otherwise explained. 



2. Methodology 

A hydrodynamic approach is used to study the CR 
effect by constructing an MHD system that incorpo- 
rates CRs as a second fluid. Cosmic rays are treated as 
a massless fluid, whose mass density is neglected, be- 
cause the mass is negligible by comparison with the en- 
ergy. The momentum spectrum of CR is neglected, sim- 
plifying the governing equations and making the com- 
putation much less intensive. The governing equations 
are, 
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in a small rectangular region of the ISM is computed 
initially in hydrostatic equilibrium with the form 

C?(z) = T(z) = To + (Thaio - ro)i[tanh(^^^) + 1], 

Z Wtr 

(6) 

where Tq is the temperature of the Galactic disk; Thah 
and Zhaio are the temperature and height of the disk-halo 
interface, and Wtr is the height of the transition between 
the disk and halo. The magnetic field is initially uni- 
formly aligned in the x direction, and independent of 
the height. Rotation of the disk is neglected. Given 
the initial temperature profile, the density, gas pressure, 
and magnetic pressure distribution at hydrostatic equi- 
librium are obtained by solving the ordinary diff'erential 
equations. These values are then input to the simula- 
tions. A velocity perturbation with a one-dimensional 
sinusoidal wave-form is applied to the gas in a few grid 
zones on the horizontal plane. The most unstable wave- 
length identified in a previous linear analysis is em- 
ployed to perturb the gas. In addition to wavelike per- 
turbations, a point perturbation that represents a single 
supernova explosion is also applied. Supernova shock 
waves are thought to be efficient accelerators of parti- 
cles with energies of up to 10^^ eV. Calculations suggest 
that at least 10% of the total supernova kinetic energy 
(10^^ ergs ) must be transformed into relativistic parti- 
cles to generate the observed Galactic CRs. 10% of the 
supernova energy is thus imparted to the CR energy. In 
the following, the ratio of plasma pressure to magnetic 
pressure is denoted a and the ratio of plasma pressure 
to CR pressure is denoted j3. 



3. Results and Conclusion 



where k = (a:|| - K_i)bb - K^dij is the diff'usion co- 
efficient tensor, p, V, Pg and jg are the mass density, 
velocity, pressure and poly tropic index of the plasma, 
respectively, B is the magnetic field, g is the external 
gravity acceleration, and Pc and jc are the pressure and 
adiabatic index of the CRs. The CR pressure Pc con- 
tributes energy to the plasma and modify the plasma 
momentum according to its gradient and the magnetic 
field. The magnetic field modifies the behavior of the 
plasma and in return aff'ect the CR via diff'usion along 
the magnetic field lines. 

The set of governing equations in Cartesian coordi- 
nates is solved. The modified Lax-Wendroff" scheme is 
applied to the MHD part and the bi-conjugate gradients 
stabilized (BICGstab) method is applied to the diffusion 
part of the CR energy equation in the same manner as 
described elsewhere [8, 9]. The temperature distribution 



Parker instability causes gas to concentrate in a valley 
of magnetic field lines and collapse toward the equato- 
rial plane. Previous investigations excluding CRs have 
demonstrated that if the external pressure is high, then 
plasma at the footpoints of the magnetic arches con- 
denses to form either clumps or filaments, or becomes 
stabilized, depending on the velocity of rotation of the 
Galactic disk. The filaments tend to align in the y direc- 
tion, perpendicular to the magnetic field line; longer fil- 
aments are formed as the strength of the magnetic field 
increases, because the magnetic pressure supports the 
piling up of gas in the vertical direction. 

Various diff'usion coefficients from k = 10 to 200 
were considered. Figure 1(a) displays the case with 
a = 1, 13 = 10 and k = 2. A comparison with the case 
without CRs (Fig. 1(b)) shows that the undular mode de- 
velops rapidly and the gas aggregates as spherical blobs. 
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Linear analysis and 2D MHD simulations by Kuwabara 
et al. (2004) [9] indicated that the growth rate of Parker 
instability declines as the coupling between the CR and 
the gas becomes stronger (meaning that the CR diffu- 
sion coefficient becomes smaller). The MHD simula- 
tions herein verify this result. 

Cosmic rays like plasma gas also accumulate near 
the footpoint of the magnetic loop (Fig. 2). For low k 
(stronger coupling between the two fluids), the CR pres- 
sure distribution is rather non-uniform, and a large CR 
pressure gradient force is exerted toward the top of the 
loop. This force impedes the falling motion of the gas 
to the bottom of the magnetic loop. Consequently, the 
growth rate of the instability is reduced. The growth rate 
normally increases with k. However, the growth is max- 
imal at /c = 32 - not k = 200. Furthermore, the change 
in the growth rate as k increases is small, because un- 
der strong difl'usion (a: > 40), the speed at which matter 
falls along the magnetic field line may exceed the speed 
of sound as revealed by a shock close to the bottom of 
the magnetic loop, and so the CR is redistributed, caus- 
ing the more uniform redistribution of the CR. There- 
fore, the contribution of CR to the gas dynamics is less 
important than in the low - CR pressure case. 

With respect to other parameters, increasing the 
transition height Wtr dampens the undular mode but 
strengthens the interchange mode. However, if a short 
wavelength is applied to excite the instability, then the 
perturbation in the x-z plane induces the interchange 
mode, but without greatly aff'ecting the undular mode. 
In the case of explosive perturbation, since the gas en- 
ergy overwhelms the CR energy, instability develops 
rapidly at first, as if strong coupling existed between the 
CRs and gas. Nonetheless, the growth in later stages 
is similar to that which occurs using wavelike perturba- 
tions. 

The morphology of the unstable gas varies consider- 
ably depending on the form and the perturbation and the 
height at which it is applied. The 3D results herein are 
similar to those in 2D cases considered elsewhere [9] 
with a perturbation on the central x-z plane. 

The simulation results herein can be compared with 
the observed properties of nearby molecular clouds, 
such as the Taurus cloud, located at 140 pc from earth. 
This cloud is a good object for examining spontaneous 
cloud formation due to Parker instability, since the 
neighborhood contains no energy sources such as mas- 
sive stars or supernova remnants that could trigger the 
collapse. Typical values of physical parameters in this 
region are evaluated. Speed of sound C = 5 x lO^cm/s 
(corresponding to temperature T = 30 K), density p = 
2 X lO^^g/cm^ (corresponding to approximately 1000 
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Figure 1: Contours of magnetic field lines for of = 1 and k\\ = 2 

particles per cubic centimeter), a scale height of H = 
CI ^jlnGp = 1.7 X 10^^ cm = 0.55 pc, time units of 
[t] = [HICs\ = 3.5 X 10^^ sec = 10^ yr, and a = 1 
are adopted, corresponding to a magnetic field strength 
of 10 yuG. Observations have revealed that the Taurus 
cloud have two to three mutually parallel cylindrical re- 
gions, which is separated by a distance of 5 pc between. 
The distance between two filaments in 3D simulations 
that exclude CRs is approximately 10 H (Chou et al. 
2000) [11], which value is consistent with a distance 
of 5 pc. A similarly dense morphology has also been 
observed in the Ophiuchus cloud. Notably, the models 
herein do not show elongated condensation since only a 
one-dimensional sinusoidal wave was used to excite the 
instability. However, the slendemess of the filaments 
may also be attributed to a weak rotation or strong CR 
eflfect, which factors will be explored in future study. 
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Figure 2: Contours of mass density, gas pressure and cosmic-ray pres- 
sure 
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